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We propose a new method to investigate the thermal properties of QCD with a small chemical potential (/i). 
The derivatives of the phase transition point with respect to ^ are computed for 2 flavors of p4- improved staggered 
fermions. We moreover comment on the complex phase of the fermion determinant in finite-density QCD. 



1. Introduction 

The study of the quark-gluon plasma in heavy- 
ion collision experiments is one of the most inter- 
esting topics in contemporary physics. To under- 
stand these experiments, precise theoretical in- 
puts of the phase transition of QCD are indis- 
pensable. Over the last several years, the numer- 
ical study of lattice QCD has been successful at 
zero quark chemical potential and high tem- 
perature In contrast, because the quark de- 
terminant is complex at /i 7^ and Monte-Carlo 
simulation is not applicable directly, study at fi- 
nite PL is still in the stage of development. 

Fortunately, the interesting regime for heavy- 
ion collisions is rather low density, e.g., fi ~ 
15MeV (/x/Tc - 0.1) for RHIC §. Therefore, an 
approach of Taylor expansion may be the most ef- 
ficient way, by computing the derivatives of phys- 
ical quantities in terms of at /i = to deter- 
mine these coefficients. Pioneering works in such 
a framework are given for free energy (quark num- 
ber susceptibility) Q and screening massQ. 

In this study, we investigate the transition 
point (/3c) at finite /i. A similar study has been 
done by l^]. In sec. ||, we propose a new method 
to compute the derivative of the physical quan- 
tities with respect to p. The result of the second 
derivative of (3c for 2 flavor QCD is given in sec. 
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H We discuss the problem of the complex phase 
in sec. ^. Section ^ presents the conclusions. 

2. Reweighting method for /^-direction 

Basically most of the attempts to study finite 
density QCD have been tried by performing sim- 
ulations at /i = and reweighting by the identity: 

(0)(^,.)=(OW^)(^„,o)/W(ft,o)> (1) 

^ ^ gaAff(lndctM(Ai)-lndotJ\/(0))g-Sg(/3)-HSg(/3o)^ 

where M is the fermion matrix, Sg is the gauge 
action, Nf is the number of flavors and a is 1 or 
1/4 for the Wilson or staggered fermion. 

The reweighting factor of the gauge part is 
easy to compute However, the fermion part 
is highly non-local and difficult to compute in 
practice. Here we consider a method which is 
applicable for small /i. We perform a Taylor 
expansion for the fermion part of the reweight- 
ing factor around /i = 0. Similarly, we expand 
fermionic observables, e.g., the chiral condensate, 
{'ipijj) = {l/V)aNi{iT:M-^). If we consider up to 
the n-th derivative of both the reweighting factor 
and the fermionic observable, we can calculate the 
n-th derivative of the physical quantity with re- 
spect to p, correctly, which can be easily checked 
by performing a Taylor expansion of each expec- 
tation value of the physical quantity. The deriva- 
tives of Indet M and trM~^ can be calculated by 
the random noise method, which enables us to 
compute on a rather large lattice in comparison 
with usual studies of finite density QCD. 
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Table 1 



JU^Pf) d^/3c/dAi^ /3c(Polyakov) d^^JdiF~ 3.650 
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0.2 



.6463(25) 
.7575(19) 



-1.41(61) 
-1.51(42) 



3.6489(27) 
3.7590(20) 



-1.71(67) 
-1-61(36) 
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Figure 1. Chiral susceptibility at to = 0.2. 



Moreover, we should notice that the odd order 
derivatives of In det M and trM~^ are pure imag- 
inary and the even order derivatives are real at 
fi — 0. This property is proved using the iden- 
tities: Mt(^) = ^5M(-/i)75, and = 
(— 1)"75^^(— /i)75. From this property, we can 



explicitly confirm that, if the operator we com- 
pute is real at fi = 0, e.g., {ipip), the first deriva- 
tive of the expectation value is zero at /i = 0, as 
we expect from the symmetry under the change 
from /i to We moreover can easily estimate 
the phase factor of the fermion determinant from 
this feature, which is discussed in sec. ||. 

3. Results for Ni = 2 improved staggered 

We calculate the second derivative of /3c with 
respect to /x. We employ a combination of the 
Symanzik improved gauge and the p4-improved 
staggered fermion actions 0. It is known that 
this action makes the discretization error of pres- 
sure small at high temperature, and Tc obtained 
by this action is consistent with that using im- 
proved Wilson fermions j|] . To include the finite 
density effect, we multiply the hopping terms pro- 
ceeding n-steps in positive and negative temporal 
directions by e"^ and e""'' respectively. 

We investigate the transition points at to = 0.1 
and 0.2. The corresponding pseudo-scalar and 
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Figure 2. Phase transition points for m = 0.1 
(upper) and 0.2 (lower). Gray lines are the fitted 
lines. 



vector meson mass ratios are mps/my ~ 0.7 and 
0.85. We compute Polyakov loop, chiral conden- 
sate, and their susceptibilities. The simulations 
are performed on a 16'^ x 4 lattice at /3 = 3.64- 
3.65 with a total of 45000 trajectories for m — 0.1, 
arid (3 = 3.74-3.80 with 76700 (Polyakov) or 69000 
(i^ip) trajectories for to = 0.2. 10 sets of Z(2) 
noise vectors are used for each trajectory to com- 
pute the reweighting factor up to the second or- 
der. We plot, in Fig.|l|, the chiral susceptibility at 
TO = 0.2. This figure shows that the peak posi- 
tion becomes lower as fj, increases. Figures^ show 
the transition point determined by the peak po- 
sition of the Polyakov loop susceptibility and the 
chiral susceptibility as a function of /i^. Because 
the first derivative is zero, as we discussed before, 
we fit the data for (3c by a straight line in /.t^, 
fixing j3c at /i = 0, in a range, fi^ < 0.005 and 
0.01 for TO = 0.1 and 0.2, respectively, in which 
the phase problem is not serious (see sec. ^). We 
then obtain the results in Table |l|. Quark mass 
dependence of d^/?c/d/i^ is not visible within the 
error. We denote the interesting value of for 
RHIC in Figs.S. The shift of /3c from ^ = is 
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Figure 3. Expectation value of cos 9 near (3c as 
a function of /i on a 16^ x 4 lattice. 



found to be small at this point. 

The second derivative of Tr can be estimated 



by 



We obtain the 



beta-function from the data of the string tension 
in Ref.§, a-\da/d(3) = -2.08(43) at (/3,m) = 
(3.65,0.1). We then find Tci<PTJd^il^^^J « -0.2 
at m = 0.1. The curvature of the transition line 
in the small fi region seems to be a little smaller 
than the phenomenological expectation, a ten- 
dency also suggested by the result of Fodor and 
Katz [|. 

4. Comment on the phase (sign) problem 

If the reweighting factor in Eq.(^ changes its 
sign frequently due to the complex phase of the 
quark determinant, the calculation of the expec- 
tation value becomes very difficult, which hinders 
the study at /i 7^ 0. The phase of the quark deter- 
minant is zero at fj, = and becomes larger as 
increases. Therefore, it is important to consider 
for which value of /i the problem becomes severe. 

As discussed in section ^, the phase can be ex- 
pressed using the odd terms of the Taylor expan- 
sion of IndetAf. Denoting detM = |detM|e*^, 



aNf Im 



d In dot M 
d/i 



For small /i, the first term, aiVflm 
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is dominant. 



Because the 



average of the phase factor is zero, the phase fluc- 
tuation is important. We investigate the standard 

deviation of -^Im tr (^/^^^^f ) • We obtain that 
the magnitude of the phase fluctuation is about 
21/i for TO = 0.1 and 16/^ for m = 0.2 near /3c. 



Consequently the phase problem starts to appear 
at ^l ^ 0.07(0.1), i.e., Uphys/Tc - 0.3(0.4) for 
TO ~ 0.1(0.2), since the phase problem arises if 
the phase fluctuation becomes larger than 0(1). 
We notice that the value of /x for which the phase 
fluctuations becomes signiflcant decreases as the 
quark mass decreases. Figure || shows the aver- 
age of the phase factor e'^. We can see that the 
average becomes small around those values of /i. 

5. Conclusions 

We propose a new method based on Taylor ex- 
pansion to investigate thermal properties of finite 
density QCD. We computed the chiral suscepti- 
bility and the Polyakov loop susceptibility for 2 
flavors of p4-improved staggered fermions, and 
found that (3c and Tc become smaller as /i in- 
creases. Our results suggest that the discrepancy 
of Tc from /i = is small in the interesting region 
for heavy-ion collisions. We also estimated the 
complex phase of the fermion determinant for a 
16^ X 4 lattice. We found that the sign problem 
is not serious in the range of ^iphys/Tc < 0.3-0.4 
for TO — 0.1-0.2, which covers the regime of the 
heavy-ion collisions at RHIC. 
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